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For systems close to equilibrium, the relaxation properties of measurable physical quantities are 
described by the linear response theory and the fluctuation-dissipation theorem (FDT). Accordingly, 
the response or the generalized susceptibility, which is a function of the unperturbed equilibrium 
system, can be related to the correlation between spontaneous fluctuations of a given conjugate 
variable. There have been several attempts to extend the FDT far from equilibrium, introducing 
new terms or using effective temperatures. Recently, Prost, Joanny, and Parrondo [Phys. Rev. Lett. 
103, 090601 (2009)] have shown that the FDT can be restored far from equilibrium by choosing 
the appropriate variables conjugate to the external perturbations. Here, we apply the generalized 
FDT to a system perturbed by time-dependent deterministic forces and acting under the influence 
of white a-stable noises. 

PACS numbers: 05.10.Gg, 05.70.Ln, 05.40.Fb. 



I. INTRODUCTION 

The fluctuation-dissipation theorem (FDT) connects 
correlation functions to linear response functions and 
constitutes a useful tool in investigations of physical 
properties of systems at thermodynamic equilibrium 
[llj l2(- By virtue of the FDT, measurable macroscopic 
physical quantities like specific heats, susceptibilities or 
compressibilities can be related to correlation functions of 
spontaneous fluctuations. For systems weakly displaced 
from equilibrium, the FDT allows one to express the lin- 
ear response of physical observablcs to time-dependent 
external fields in terms of time-dependent correlation 
functions. Accordingly, departures from the FDT can 
be expected for far-from equilibrium situations and have 
been demonstrated in various aging, glassy and biological 
media 

On the other hand, the wealth of theoretical, experi- 
mental and numerical research indicate that the FDT is 
a special case of more general fluctuation relations that 
remain valid also in a specific class of nonequilibrium sys- 
tems Following a former generalization of FDT 
15| based on the identity derived by Hatano and Sasa 
1 61 ] , we discuss here an extension of the fluctuation the- 
orem to stochastic models obeying Markovian dynamics 
and driven by white a-stable noises. We apply the gener- 
alized fluctuation-response theorem to this case and ana- 
lyze the regime in which linear response theory becomes 
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invalidated. We illustrate our results with the simple ex- 
ample of an oscillator coupled to a non-equilibrium bath 
whose action is represented by a Levy white noise. 

Let us review the FDT introduced in [ll| and extended 
to arbitrary observables in Q. The theorem applies to 
any Markov process x(t) whose dynamics depends on a 
set of parameters A. We study the linear response of 
the system to perturbations \(t) = \q + 5X(t) around a 
reference state, p ss (x;Xo) which is the stationary state 
corresponding to constant parameters Ao- Given an ar- 
bitrary observable A(x), the response can be written as 

(A(t)) - (A) ~ f XA n {t - t')5\ 7 (t')dt', (1) 



where A{t) = A(x(t)) and the brackets (. . . )o indicate 
an average over the reference state p ss (x; Ao); summation 
over the repeated index 7 is assumed and XA,-y is the 
time-dependent susceptibility of variable A with respect 
to variations of A 7 (i.e. perturbations of the 7 component 
of A). The FDT relates this susceptibility to correlations 
measured in the reference unperturbed state 



XA, 7 (t-t') = -(A(t)X 1 (t')) , 



(2) 



where X^(x) is the variable conjugate to the perturbation 
A^ which is defined as 



X 7 {x) 



din p ss (x; A) 



5A„ 



(3) 



A=A 



In this definition p ss (x;\) is the stationary probability 
distribution function (PDF) for constant values of the 
parameters A. 
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If the reference state is the Gibbs equilibrium state 
corresponding to a temperature kT = P^ 1 and a Hamil- 
tonian H(x; A), the stationary PDF p BB (x; A) assumes the 
form p BS (x;X) = exp[—/37i(x;X)]/Z(/3,X) and the conju- 
gate variable reads 



XJx) 



1 d H(x; A) - F(P, A) 



kT 



dXr, 



(4) 



A=A 



where F = —kThiZ stands for the free energy. Accord- 
ingly, A" 7 can be interpreted as the fluctuation of the 

dH(x;\) 



quantity dH ^ o) 



A=A 



XJx) 



1 

kT 



d'H{x;Xo) / dH(x;X ) 



dX~, 



dX-, 



(5) 



For instance, if A 7 is a force coupled to a coordinate x 7 , 
i.e., if the control parameter appears in the Hamiltonian 



as — A 7 x 7 , then the conjugate variable AT 7 



(x 7 ))/(kT) represents fluctuations of x 1 . 

On the other hand, if the reference state p ss (x;Ao) is 
not an equilibrium state, the conjugate variables defined 
by Eq. ^ do not have any straightforward physical inter- 
pretation [H[l5|. In this Letter, we apply the generalized 
FDT to a system obeying non-equilibrium Markovian dy- 
namics driven by Levy white noise. The system of this 
type may be conceived as a generalization of Brownian 
motion: the particle undergoing Levy superdiffusion is 
performing motion with random jumps and step lengths 
following a power-law distribution. As a result, the width 
of the distribution of particles grows superlinearly with 
time [13, [H| signaling anomalous dynamics. A charac- 
teristic feature of Levy flights is the dominance of rare 
but large jumps. This type of anomalous transport has 
been found ubiquitous in nature [13, [HI and serves as a 
suitable model for various physical phenomena like atmo- 
spheric turbulence [20| . transport in turbulent plasmas 
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activation kinetics in non-equilibrium baths [22j], 
transport in fractured materials 1231 . epidemic spread- 
ing [24[ , dispersal of banknotes 12511 or light scattering in 
heterogeneous dielectric media [26| . 



II. LINEAR SYSTEM DRIVEN BY LEVY 
WHITE NOISES 

We proceed to discuss response properties of an over- 
damped Levy-Brownian particle moving in a parabolic 
potential that is subject to a deterministic time- 
dependent force field f(t) and a white Levy noise 
resulting from the fluctuating environment. The corre- 
sponding Langevin equation reads 



x(t) = -ax + f(t) + C(t) 
x(0) = x 



(0) 



The white Levy noise ((t) is defined as the time derivative 
of a stationary Levy process [I3,[13, i- e -, the integral over 
time 



({s)ds = z(t) 



(7) 



represents a stochastic process with independent incre- 
ments whose probability density p a ^(z,t) is a stable 
Levy distribution. Consequently, the Fourier trans- 
form of the probability density (characteristic function) 
f(k,t) = e lkz ^p a ^{z,t;a ,iJ,)dz reads 



tp(k, t) = exp ikp t — <Tq \k\ a t (l — if3 sign(fc) tan 
for a / 1 and 



7TQ!\ 

~2J 



(8) 



ip(k, t) = exp 



ikji^t — uo\k\t 1 



if}— sign(fc) In \k 

7T 



(9) 

for a = 1 [28|. Here a £ (0, 2] is the Levy (stability) in- 
dex, /3 £ [— 1, 1] is the skewness parameter (for j3 = the 
distributions are symmetric), uq > represents the noise 
intensity, and pg £ K is a location (shift) parameter. The 
Gaussian distribution is a special case of symmetric sta- 
ble distribution with a = 2 and (3 = 0. In this case pot is 
the mean of the Gaussian random variable z(t) and a^t 
is its variance. For a < 2, stable probability densities 
exhibit long tails and divergent moments: the asymp- 
totic (large z) behavior of the corresponding PDF is then 
characterized by a power-law p a ,p(z, t; uq, p) cx |z|^' 1+Q \ 
Under those circumstances, equation ([6]) is associated 
with the space-fractional Fokker-Planck-Smoluchowski 
equation (FFPE) JH-iH: 



dp(x, t) 



dt 



d 

— [p - ax + f(t)]p(x,t) 



(10) 



+ag-Qr-r-p(x,t)+a%pt<m 



wa d d 



a-l 



2 dx d\x 



Here, the fractional (Riesz-Weyl) derivative is defined 

by its Fourier transform F j^ a ^ a f(x) = — \k\ a F [f{x)\ 

HHH. Accordingly, Eq. $TU$ has the following Fourier 
representation 

= ~ ak §k m t] + lh [Mo + /(t)] m t] 

7TCX 1 

1 — i(3 sign(fc) tan — p(k,t) 7 

where p(k, t) = F [p{x, t)]. In what follows, we adhere to 
the analysis of strictly a-stable random variables [27j, i.e. 
ones for which p = 0, and additionally /3 = if a = 1. 

Since our original Langevin equation Q is linear, its 
solution depends linearly on the stable process L a ^{t). 
Accordingly, the probability density of the solution, 
p(x, t|xo, 0), has the form of an (a, /3)-stable Levy distri- 
bution with time-dependent location p(t) and scale cr(t) 
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parameters [33[. By analogy, its characteristic function 
is given by (cf. Eqs. © and ©) 



p(k,t) = cxp ikfi(t) - <r a (t)\k\ a (l - i/3sign(fc) 



tan — 
2 / 



We can now insert this ansatz into the Fokker-Planck- 
Smoluchowski equation (JTTJ. Since the derivative with 
respect to k appears multiplied by k in the non 

analyticity of \k\ a at k = does not create any singularity 
in the equation. The real part of Eq. (|TTj) yields the 
following evolution equation for the scale parameter o~{t) 



whereas the imaginary part gives 
[fi + afi- f(t)]k= [-a<7 a ~ ] 



(12) 



(13) 



x /3tan — |fc| Q sign(fc). 



The right hand side of Eq. (fTB"]) vanishes due to Eq. 
From the left hand side one gets the evolution equation 
for the location parameter p,(t): 



= -a/i + / (f). 



(14) 



The evolution equations (|T2"j) and (IT41 are completed with 
the initial conditions fi(0) — xq and <r(0) = (we are 
calculating probability densities conditioned to x(0) = 
xq). The solution of these differential equations are 



/i(t) = 



x 



e as f{s)ds 



and 



<y(t) = a 



— (1 - e - aat ) 
aa 



l/a 



(15) 



(16) 



where <7o is the scale parameter of the correspond- 
ing a-stable density. For a constant force f(t) = f, 
the long time asymptotics of the above equations are 
lim^oo n(t) = f/a and lim^oo a(t) = a /(aa) 1/a . 



III. THE CONJUGATE VARIABLE 

To determine the conjugate variable to the external 
force, we need the stationary distribution in position 
space for a constant force /. Despite the characteristic 
functions of stable distributions assume closed expres- 
sions, the corresponding PDFs have a known simple an- 
alytical form [23, HH only in a few cases: For a = 2 
and /3 = the resulting distribution is Gaussian; for 
a = 1, /3 = one gets the Cauchy distribution; finally, 
for a = 1/2, (3 = 1 one gets the Levy- Smirnoff distribu- 
tion. In this section, we give explicit expressions for the 
conjugate variable for these three cases. 



For a = 2 and (5 = 0, the time dependent solution of 
the corresponding Langevin Eq. © is 



P2,o(x,t\x ,0) 



1 



exp 



2a 2 (t) 



(17) 



with (i(t) and cr(t) given by Eqs. (|T3J) and ([15]). The 
stationary solution p ss (x) for a constant force / is ob- 
tained by replacing p,(t) and cr 2 (t) by their stationary 
values, f/a and <7g/(2a), respectively. We then get the 
non-equilibrium potential 



4> = - h\p ss (x) 



1 



In 



a a(x — f /a) 



(18) 



""0 u o 
and the conjugate variable can be easily derived as 

9f /=0 



2.T 



(19) 



which is proportional to expected, since the Gaus- 
sian case corresponds to a Brownian particle in equilib- 
rium. 

For a = 1 and j3 = 0, the time dependent solution of 
the corresponding Langevin Eq. ([6]) is the Cauchy distri- 
bution 



Pi fi (x,t\xo,Q) 



a(t) 



1 



7T [aj-/i(t)r + o- 2 (t) 



(20) 



and the stationary solution for a constant force / is ob- 
tained replacing p,(t) and a(t) by their stationary values, 
f/a and o~o/a, respectively. The non-equilibrium poten- 
tial in this case is given by 



= - In 



+ \n[(x-f/a) 2 + (a /a) 2 



so that the conjugate variable takes the form 

2x 

An 



a [a: 2 + (ao/a) 5 



(21) 



(22) 



which is proportional to x only for small values of x and 
becomes proportional to 1/x for large x. This large x 
behavior ensures the convergence of all the moments of 
Xc, whereas for the Cauchy case \x\ v exists only if v < 1, 
see [13 

Finally, for a = 1/2 and /3 = 1 the solution to Eq. © 
is the Levy-Smirnoff PDF 



P 1 /2,l(x,t\x ,t Q ) 



a(t) 



2ir [x — fi(t)} 



■ exp 



a(t) 



2{x-p{t)) 
(23) 

for x > p,(t) and Pi/2,i(x,t\xo,0) = for x < The 
stationary values of /z(t) and cr(i) are in this case f/a and 
4<7o/a 2 , respectively. Inserting these values, one can eas- 
ily obtain the non-equilibrium potential <fi = — lnp ss (x) 
and the conjugate variable 



4(To — 3a x 
2a 3 x 2 ' 



x > 0. 



(24) 
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IV. SUSCEPTIBILITY AND RESPONSE 

The main objective of the current work is to compare 
the response of the system to external perturbation as 
calculated directly from the definition 



(*(*)> 



X(x)p(x, t)dx, 



(25) 



or, otherwise determined by the generalized susceptibility 
x(t) = ^(X(t)X (0))o within linear response theory: 



LR 



X (t - s)f(s)ds. 



(26) 



We restrict our analysis to the Cauchy case, a = 1, /? = 
0, and denote the conjugate variable as X = Xc, with 
Xc given by Eq. (|2"2")l . In this case, the time-dependent 
average (|25|) can be calculated exactly with the probabil- 
ity density: 



p(x,t) 



p{x,t\x o ,0)p(x o )dx o 



where 



O"0 



an xr. 



(27) 



(28) 



and p(x, t\xo, 0) is given by Eq. (|20j). 

On the other hand, the FDT relates the susceptibility 
with the autocorrelation of the conjugate variables in the 
reference state, i.e., for / = 0. The autocorrelation is 
defined as 



(X(t)X(Q)) = 



2.r 



2y 



a[x 2 + (aa/a) 2 } a[y 2 + (ao /a) 



n[(x-Kt)) 2 + v 2 (t)] 

TO 

air [y 2 + {<r /a) 2 



■dxdy 



where fi(t) = e~ at y and a(t) = cr [(1 - e~ at )/a]. The 
final result is surprisingly simple: 



(X(t)X(0)) 



1 

2$ 



(29) 



From the above, the generalized susceptibility can be de- 
rived by differentiation with respect to time (see Eq. ©): 



x ( t ) = ±(X(t)X(0)) = --^ 



dt 



2a 2 



(30) 



In further calculations, for the sake of simplicity, it is 
assumed that a = 1 and <To = 1, so that (X(t)X(0)) = 



and x(t) 



■ie-* 



In order to test the linear response theory for our dy- 
namic Markov system subjected to Levy white noise, we 
calculate the response of the conjugate variable X to two 
different time dependent perturbations: the sum of a 



A 
x 
V 



0.2 



-0.5 



Irt 

exact 




150 200 




200 



FIG. 1. Response of (X(t)) to external drivings fi(t) = 
sin(t)/10 + f/100 (upper plot) and f 2 (t) = tsin(t)/100 (lower 
plot). The solid and dotted lines present an exact result 
(Eq. (|25|l ) and a result constructed by use of the linear re- 
sponse theory (Eq. (|26|l 'l. respectively. 



small periodic and a linearly increasing force, fi(t) — 
sin(t)/10 + i/100; and a periodic force with increasing 
amplitude, fo(t) = isin(t)/100. Figure [1] displays the 
exact evolution of (X(i)) and the result obtained from 
the linear response theory. For small perturbations, (i.e. 
short times) in both cases, linear response theory yields 
an accurate estimation of the response. In the case of 
fi(t), for large times the response (X(t)) is insensitive 
to the sinusoidal component of the force, which is small 
compared with the linear part. This is due to the pecu- 
liar form of the conjugate variable X given by Eq. (|22[) . 
For a constant force /, the mean value of X is 



(X) 



<T0 

an 



dx 



2.r 



[x - f/af + K/a) 2 a[x 2 + [a a /af 



2/ 



P + 



Ao 2 



which yields (X) = 
t = 200, f x (t) ~ 2). 



-0.5 for cr = I and / = 2 (at 
A similar saturation effect is not 
observed for the sinusoidal force fait). 

We can apply the FDR to any function A(x) with fi- 
nite average. Due to the long tails of stable distributions, 
only moments with v < a converge (v < 1 in the 
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case of Cauchy distributions) [27|. Moreover, those mo- 
ments are even functions of x and, for symmetry rea- 
sons, the correlation with X vanishes: (\x(t)\ v X (Q)) = 0. 
Consequently, the deviation of (\x(t)\") with respect to 
its reference value is non linear in the perturbation f(t). 
On the other hand, we can obtain non-trivial results for 
odd fractional moments A(x) = sign [x(t)] [x^t)^, whose 
average in the reference state vanishes (^4)o = 0. The 
corresponding correlation function reads 



(A(x(t))X(0)) 



sin [ttv /2] 



and the generalized susceptibility is given by 

XA{t) 



sin [itv/2] 



(31) 



(32) 



In the spirit of the former definition, see Eq. ([25)1 . the 
exact value of (A(x(t))) can be calculated as 



(A(x(t))) = 



A(x(t))p(x,t)dx (33) 
A(x(t))p(x, t\xo, 0)p(xo)dx, 



where p(xo) and p(x,t\xo,0) are given by Eqs. (j2TJ)) 
and l[2"8|) respectively. Figure displays the compari- 
son of the exact evolution (sign [x(t)] ^(i)! 1 / 2 ) with the 
linear response approximation J^XA^t — s)f(s)ds. As 
previously, linear response theory is valid for weak per- 
turbation up to / ~ 0.5. However, for the fractional 
moment and the linearly increasing force (upper plot), 
we do not observe saturation. 





200 



200 



FIG. 2. Response of the (sign [x(t)] \x{t)\ 1/2 ) to f(t) = 
sin(i)/10 + 4/100 (upper plot) and f(t) = tsin(t)/100 (lower 
plot). The solid line presents exact results, see Eq. (|33p . while 
the dotted line results constructed by use of the linear re- 
sponse theory, see Eqs. (|26jl and pip . 



V. SUMMARY AND CONCLUSIONS 

We have shown that the generalized FDT can be ap- 
plied to linear systems driven by Levy noise. The FDT 
allows one to calculate the susceptibility of any observ- 
able and then the response to any small time dependent 
perturbation. For a noise distributed according to the 
Cauchy distribution, we have calculated the susceptibil- 
ity of the conjugate variable Xq and the susceptibility of 
odd fractional moments (sign[x(t)}\x(t) 1 ' |), which have a 
simple exponential behavior. From these susceptibilities 
it is easy to get simple analytical expressions for the re- 
sponse of the system using Eq. ([26]) . We have to notice 
that, although the exact response can be calculated an- 
alytically using Eq. (|2"5"|) . the corresponding integrals are 
cumbersome and can be only solved numerically in the 
simplest cases. Therefore, the generalized FDT is shown 
to be a useful analytical tool to deal with these type of 
systems. 

It is still not obvious whether the conjugate vari- 
ables that we have calculated for the Cauchy and Levy- 
Smirnoff noises have any physical meaning, besides the 
one provided by the generalized FDT itself. The general- 
ized FDT shows that these conjugate variables represent 



the change in the probability distribution of the system 
under the perturbation. In equilibrium, this change is 
also related with the energy that the system absorbs from 
the perturbation. On the other hand, for non-equilibrium 
systems, the lack of conserved quantities prevents such an 
interpretation. For instance, in the- case of the- harmonic 
oscillator driven by a Cauchy-Levy noise, (x) and higher 
moments diverge (cf. (35{). Consequently, both the po- 
tential energy of the system in the harmonic potential 
and the work done by the external force f(t) also diverge. 
The system is plagued by divergent quantities. However, 
the conjugate variable Xq given by Eq. (f2"2"j) has finite 
moments and still captures the dynamical response of the 
Levy particle. Summarizing, although systems driven 
by a-stable noises might significantly differ from their 
Brownian (equilibrium) counterparts j36M38| due to their 
heavy tail asymptotics, we have shown that in such far- 
from-equilibrium situations some concepts from weakly 
perturbed equilibrium systems can be still used. 

One of the drawbacks of the generalized FDT derived 
in fl5| is the difficulty to find the conjugate variable, since 
it requires the knowledge of the stationary state. We have 
been able to find this stationary state for a linear system. 
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An interesting open question is whether this state, or 
some slight modification, can still be used to calculate 
susceptibilities in the presence of weak non-linearities. 
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